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ABSTRACT 

Aims. The electron plasma frequency oj pe and electron gyrofrequency Q e are two parameters that allow us to describe the properties of a plasma 
and to constrain the physical phenomena at play, for instance, whether a maser instability develops. In this paper, we aim to show that the maser 
instability can exist in the solar corona. 

Methods. We perform an in-depth analysis of the ci> pe /Cl e ratio for simple theoretical and complex solar magnetic field configurations. Using the 
combination of force-free models for the magnetic field and hydrostatic models for the plasma properties, we determine the ratio of the plasma 
frequency to the gyrofrequency for electrons. For the sake of comparison, we compute the ratio for bipolar magnetic fields containing a twisted 
flux bundle, and for four different observed active regions. We also study how co pe /Q. e is affected by the potential and non-linear force-free field 
models. 

Results. We demonstrate that the ratio of the plasma frequency to the gyrofrequency for electrons can be estimated by this novel method 
combining magnetic field extrapolation techniques and hydrodynamic models. Even if statistically not significant, values of (L> pe /Cl e < 1 are 
present in all examples, and are located in the low corona near to photosphere below one pressure scale-height and/or in the vicinity of twisted 
flux bundles. The values of aj pe /Cl e are lower for non-linear force-free fields than potential fields, thus increasing the possibility of maser 
instability in the corona. 

Conclusions. From this new approach for estimating co pe /£l e , we conclude that the electron maser instability can exist in the solar corona above 
active regions. The importance of the maser instability in coronal active regions depends on the complexity and topology of the magnetic field 
configurations. 
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1. Introduction 


Loss-cone driven instabilities play an important role in space 
plasmas. The emission of non-thermal radiation is still a puz¬ 
zling process despite the accurate observations made in other 
fields, such as planetary magnetospheres dTreumannj 1120061) . 
This emission is often explained by the cyclotron maser mech¬ 
anism for which the X-mode and its harmonics provide the 
escaping radiation. The growth rate of the X-modes depends 
on the value of the ratio of the electron plasma frequency to 
the ele ctron gyrofrequency (hereafter E e ). ISharma & Vlahosl 
(1984) studied the development of the X-modes and O-modes 
under the physical conditions of plasmas in the low corona. 
The study has been performed for values of E e less than 2.5. 
The authors showed that (i) the electron maser instability dom¬ 
inates the emission for E e < 1 over other loss-cone driven in¬ 
stabilities, and (ii) the first harmonic of the X-mode dominates 
for E e < 0.35. In solar physics, the electron maser instability 
has been studied to explain the spiky radio burst observed dur¬ 


ing flares ( Holman et aPll980 ; Melrose & Dulk|[l982h . Similar 


0.35 and 1 for the growth of X-modes (e.g. Melrose et al. 

1984; Winglee 

1985; Vocks & Mann 2004; Tang & Wu 2009; 

Lee et al. 2013 

). The question arising from these studies is, Can 


electron cyclotron maser emission be a viable mechanism in the 
solar corona ? To tackle this question, we estimate the ratio E e 
in the low corona by assuming that the corona is in equilib¬ 
rium and thus can be described by a force-free equilibrium for 
the magnetic field and a hydrostatic equilibrium for the plasma 
parameters. 


The electron population of a plasma can be described 
by two characteristic frequencies amongst other parameters: 
the electron plasma frequency, oj pe , and the electron gyrofre¬ 
quency, Cl e . The plasma frequency of a cold plasma (neglecting 
the thermal velocity) is given by 
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Fig. 1. Spatial distribution of E e < 1 for the bipolar field, (top row) top view; (bottom row) side view for the potential (left) and 
niff (right) magnetic fields with varying gravity. The values of E e are averaged along the 3rd dimension. Only the E e < 1 values 
are colour-coded. 


where n e is the electron number density, e is the electric charge 
(1.6022xl0 -19 C), m e is the electron mass (9.1094xl0 -31 kg), 
and £o is the permittivity of free space (8.8542xl0 -12 F m -1 ). 
The gyrofrequency of electrons is given by 

eB 

O, = —, (2) 

m e 

where B is the magnetic field strength. The ratio of the two 
frequencies, E e = aj pe /£l e , can thus be expressed as the ratio of 
the speed of light to the Alfven speed, 


1 c 
1.147? vk’ 


(3) 


where R is the mass ratio between protons and electrons (7? = 
1836) and the coefficient 1.14 is obtained for a coronal plasma 
with a mean atomic weight of jl = 0.6. Thus, knowing the 
Alfven speed from the low-order magnetohydro static model 
developed bv [Regnier et al.l (2 008k we derive E e for different 


magnetic field configurations. This new method is then used to 
determine whether the maser instability is a viable mechanism 
in the solar corona. 

The paper is organised as follows: after summarising the 
magnetic field and hydrostatic models in Section [2j we present 
the estimate of aj pe /Q e for a simple bipolar configuration in 
Section [3j and for four different active regions in Section [4] We 
summarise our findings and conclude in Section 0 


2. Magnetic field and density models for the solar 
corona 


As in [Regnier et al. (2008 ). we combine two models to describe 
the properties of magnetised plasma in the solar corona: the 
coronal magnetic field is assumed to be a force-free field, and 
the plasma properties are derived from a stratified atmosphere 
in hydrostatic equilibrium. In the following sections, we briefly 
summarise the different numerical models. 
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2.1. Potential and non-linear force-free models 


Magnetic field extrapolations are well-known techniques for 
describing the 3D nature of the coronal magnetic fields (see 
review bv [Regnier 2013 ). Under coronal conditions, the mag¬ 
netic forces dominate the pressure gradient and gravity, and 
so we regard the corona above active regions as being well 
described by the force-free approximation (see e.g. recent re¬ 
views bv [Wiegelmann & Sakuraill2012l: lRegnieJ2Q 13l and ref¬ 
erences therein). Throughout this article, the coronal magnetic 
configurations are computed from the non-linear force-free 
(niff) approximation based on a vector potential Grad-Rubin 
( 1958 ) method by using the XTRAPOL code dAmari et al 


1997L119991) . The niff field is governed by the following equa¬ 


tions, 


V x B = aB, 


B • Vcr = 0, 
V B = 0, 


(4) 

(5) 

( 6 ) 


where B is the magnetic field vector in the domain Q above the 
photosphere, <5Q; a is a function of space defined in Cartesian 
coordinates as the ratio of the vertical current density, J z \ and 
the vertical magnetic field component, B z \ 


1 Id B y d B x \ 
B z \dx By )' 


(7) 


From Eq. (0, a is constant along a field line, but varies across 
field lines. The full descriptio n of t he Grad-Rubin iterative 
scheme can be found in 1 Amari et aL ( 199 71), and the vector 
magnetograms used as boundary conditions for the follow 


ing examples have been detailed, for instance, in lRegnier et al 

120081 ) . 

For the sake of comparison, we compute both the po¬ 
tential and niff fields following the technique developed by 


Grad & Rubin dl958l) . 


where H = kffPKprn p go) is the pressure scale-height (ks = 
1.38 10 -23 J K _1 , p = 0.6 for a fully ionised coronal plasma, 
m p = 1.67 10 -27 kg, and go = g(R©) = 274 m s -2 ), and po and 
po are characteristic values of the plasma pressure and density 
at zo- The pressure scale-height and the density po are the two 
free parameters of the model. Typical values are H = 50 Mm 
giving a coronal temperature of 1 MK and the number density 
correpsonding to po is no = 10 9 cm -3 . 

For a gravitational field varying with height, the Alfven 
speed is thus given by 


v A (x,y,z) 


g^expf ^ 

VPoPo \2H ( R 0 + zo) \R 0 + z)) 


( 10 ) 


where B is the magnetic field strength. 



Fig. 2. Distribution of E e for the bipolar magnetic field within 
the entire computational volume (cutoff at E e = 30) for the po¬ 
tential (dashed line) and niff (solid line) magnetic field models. 
The vertical dashed black line is the E e = 1 threshold. 


2.2. Hydrostatic model 

In order to define the thermodynamics parameters of the coro¬ 
nal plasma, we assume that the corona is an isothermal atmo¬ 
sphere satisfying the hydrostatic equilibrium 

- Vp +pg = 0, (8) 


In summary, for all cases (bipolar field, and active regions), 
we use two different models (potential field and niff field both 
with a varying gravity) to analyse the variation of the (L> pe /£l e 
ratio. 

3. Bipolar fields 


where p is the plasma pressure, p is the density, and g 
is the gravitational force. In agreement with the Harvard- 
Smithonian model of the solar atmosphere, we consider the 
hydrostatic equilibrium to be a reasonable assumption above 
the photosphere at a height (zo) of about 5 Mm. As shown 
in lRegnier et al.l (2 0081) . it is more appropriate to consider the 
variation with height of the gravitational field in order to satisfy 
the continuity with the properties of the solar wind. Therefore, 
we will only consider the case of varying gravity satisfying the 
following equations for the plasma pressure and density, 


| p(z ) 
I p(z) 


= Po exp 
= Po exp 


/ 4 | 

f z-zq n 

\ H (R e +z o) 1 

Uo +z)J’ 

( R ° 


l H(R q +z 0 ) 

Uo +zJJ 


and 

for g = g (z)=^, 


(9) 


As a test case for the above method, we use a simple bipolar 
magnetic field distribution on the photosphere. For all models, 
the vertical component of the magnetic field, B z , is given by 
a Gaussian distribution. To compute the niff field, the vertical 
electric current density is injected in the magnetic configuration 
following a photospheric distribution of the form 

J z (r) = 2 J z0 [r 2 - C 0 ] exp(11) 


where r is measured from the centre of the magnetic polarity 
and Co is a constant that ensures a zero net current. This J z dis¬ 
tribution is typically a second-order Hermite polynomial func¬ 
tion that allows for return current on the surface of the flux tube 
(see lRegnIe3l2012l) . The vertical magnetic field strength B z is 
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Fig. 3. Same as Fig.[T]for AR8151 


2000 G at the centre of the polarities and J z has a maximum 
strength of 10 mA m~ 2 . Both the potential and niff models are 
computed on a 148x148x148 computational box with a pixel 
size of 1 Mm. 

In Fig. E we plot the spatial distribution of E e < 1 for 
the potential and niff fields with varying gravity (from left to 
right). The first row shows a top view (xy-plane), whilst the 
second row shows a side view (xz-plane). The values of E e are 
averaged along the third dimension. The values of E e < 1 are 
localised where the magnetic field strength is strong, below a 
height of 25 Mm. For this bipolar magnetic field configura¬ 
tion, the values of E e are similar across both models (see also 

Tabled]). 

The distributions of E e within the volume are shown for 
both models in Fig. [2| (values of E e between 0 and 30). The two 
distributions have similar shapes with two main peaks around 
10 and 15, and a tail for large values of E e which tends to zero. 
For E e < 1, the distributions have a number of occurrences 
which is not significant whatever the magnetic field model. We 
confirm this in Tabled] by listing the percentage of E e in given 
interva ls. We use the same intervals as in Isharma & Ylahosl 
(1984). The percentage of E e < 1 is about 1.5% for all models, 
whilst ~95% of E e values are above 2.5. In this example, the 
models with varying gravity give similar results whether or not 
the magnetic configuration contains electric currents. Thus, the 
maser emission is weakly influenced by the electric currents or 



0 50 100 150 200 250 300 

x (Mm) 



the twist and shear in a simple magnetic configuration (without 
a topology). 

The conclusions from this test case are (i) values of E e < 1 
exist in coronal-like magnetic field configurations, and (ii) 
these values are located above large photospheric magnetic 
field strengths. 


4. Active regions 

We now analyse the values of co pe /£l e for four active regions at 
different stages of their evolution. 

AR8151 

The active region AR8151 is an old decaying active region with 
moderate and highly twisted flux bundles observed as a fila¬ 
ment and an X-ray sigm oid, respectively dRegnier et alJl2002t 
Regnier & Amaril2004l) . In Fig. [3] the E e < 1 values are located 
where the magnetic field strength is large, i.e. above sunspots 
and plages. For the potential model, the values of E e < 1 
are concentrated below 30 Mm, i.e. within one pressure scale- 
height above the surface. For the niff model, the spatial distri¬ 
bution shows significant values of E e above 50 Mm; these val¬ 
ues characterise the influence of electric currents present in the 
magnetic configuration and are associated with highly twisted 
flux bundles and sheared arcades. The twisted flux bundles ex- 
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cj„ ,/Q, 

Fig. 4. Same as Fig.|5]for (a) AR8151, (b) AR8210, (c) AR9077 

isting in this active region generate a local increase in the mag¬ 
netic field strength, i.e. in the local Alfven speed (see Eq.fTOl). 
As can be seen in Fig. [4k, the distribution of E e peaks at 2 
for the niff model and at 5 for the potential model. Both distri¬ 
butions are narrower (with a typical width of about 2) than in 
the magnetic bipole example presented in Section[3] Unlike the 
bipolar field, there is a clear difference between the potential 
and niff models; the E e values are smaller for the niff magnetic 
field. Both distributions have a tail that becomes statistically 
insignificant for values of aj pe /Qe above 10. These results lead 
us to consider that, for this active region, the niff model has a 
different behaviour than the potential field model owing to the 
presence of twisted flux bundles and the complexity of the mag¬ 
netic field; a proxy of the complexity of the coronal magnetic 
field is given by the complexity of the photospheric magnetic 
field. In terms of statistics, TableQ] shows that values of E e < 1 
represent -0.8% for the potential model and -1.4% for the niff 
model. The main difference is in the percentage of E e between 
1 and 2.5, which reaches 40% for the niff field, whilst it is just 
around 8% for the potential model; as noticed previously, this 
is due to the strong electric currents present in the active region, 
which imply sheared and twisted magnetic field lines leading 
to the conclusion that more realistic magnetic fields like the niff 
model are more favourable to maser emission. 



w pe /Q e 



0 10 20 30 

u p ,/0. 


and (d) AR10486 
AR8210 

AR8210 is a relatively new active region which has pro¬ 
duced numerous flares. The main features in AR8210 are a 
clockwise-rotating sunspot surrounded by diffuse polarities of 
opposite sign, and a small emerging parasitic polarity inter¬ 
acting with the pre-exis ting magnetic topology of the field 
jRegnier & Canfiel32006i) . A time series of niff magnetic field 
extrapolations has shown that the 3D magnetic configuration 
is close to potential (no twisted bundles, and small noticeable 
sheared arcades) and does not contain a large amount of free 
magnetic energy (estimated to A E m = 2.6 10 31 erg), i.e. not 
sufficient to produce large flares. 

In Fig. 0 we plot the spatial distribution of E e within the 
whole computational volume. As for the bipolar field described 
in Section [3l the values of E e < 1 are located just above the 
strong magnetic field regions and do not extend above 30 Mm. 
In Fig. the distributions of E e indicate that the significant 
values of the ratio are above 3. As noted for AR8151, the niff 
model decreases significantly the values of E e compared to the 
potential field model. The E e distribution peaks between 7 and 
13 for the niff model, whilst it peaks at 27 for the potential 
model. In terms of statistics (see Table []]), the percentage of 
E e < 1 is less than 0.5% for both models, and the percent¬ 
age of H, > 2.5 is about 98%. These results are consistent 
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Fig. 5. Same as Fig.[T]for AR8210 


with the fact that no twisted flux bundles or highly sheared ar¬ 
cades have been found in the 3D extrapolated magnetic field 
dRegnier & Canfiel3l2006l) . 


AR9077 

AR9077 produced a X5.7 flare on 2000 July 14 known as 
the Bastille Day flare (e.g. iKosovichev & Zharkoval 2001 


Wang et alJl2005h . The magnetic field extrapolation was per¬ 


formed after the occurrence of the X-class flare and thus de¬ 
scribes the relaxation phase associated with post-flare loops. 
It has been shown that the magnetic field configuration is quite 
close to a potential field with a small amount of shear and twist. 
These magnetic field properties are reflected in the spatial dis¬ 
tribution of E e (Fig. ft: the E e < 1 values are concentrated 
above regions of large magnetic field strengths as was seen for 
the previous examples. Both potential and niff fields have simi¬ 
lar spatial distributions below 30 Mm. Futhermore, the niff field 
exhibits an excess of low E e values above 30 Mm indicating the 
presence of twisted/sheared magnetic field lines increasing lo¬ 
cally the magnetic field strength. Compared to AR8151 (see 
Fig. [3]), the volume containing the E e < 1 values above 30 Mm 
is small, and does indicate a small amount of twist or shear 
in the magnetic configuration. In Fig. [4]c, the distributions for 
both models peak at about 1 with a percentage of values with 


E e < 1 between 14% and 19% (see Table []])• It reinforces the 
fact that the magnetic configuration of AR9077 at this stage of 
its evolution is close to potential, although the number of lo¬ 
cations where the maser instability could occur is large (about 
20%), filling the volume below 30 Mm almost entirely. 


AR10486 

The active region AR10486 is the source of the series of strong 
flares in October-November 2003 known as the Halloween 
events, including the strongest flare recorded for the corre¬ 
sponding solar cycle estimated to be a X28 flare on November 
4. The magnetic field configuration presented here is obtained 
before the occurrence of the XI7 flare located near the disk cen¬ 
tre on October 28. For this flare, the energetics and magnetic 


topology have been studied in detail (see e.g. iRegnieretal 
l2005k Mandrini et al.ll2006t IZuccarello et al.l [2009h : the mag¬ 
netic energy contained in this active region has been estimated 
to be above 10 33 erg, consistent with a large total photospheric 
magnetic flux. For the potential model (see Fig. O, the E e < 1 
values are located above strong photospheric magnetic field 
strength and extend in a volume below 70 Mm. For the niff 
field, the model generates low E e values high in the corona up 
to 170 Mm. These high-altitude, low values of E e correspond 
to several pressure scale-heights. This indicates that not only is 
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Fig. 6. Same as Fig.[l]for AR9077 



the pressure scale-height important for this model (as noted for 
the previous active regions), but also that the magnetic scale- 
height is relevant when the total unsigned magnetic flux on the 
photosphere is large. In Fig.|4]i, all the distributions of E e peak 
at 1; however, the niff distribution has a larger number of oc¬ 
currences and a narrower width, with 80% of the values be¬ 
ing between 0 and 2.5, compared to the potential distribution 
which has an extended tail towards large values of E e . This is 
reflected in Table [TJ which shows that the niff magnetic con¬ 
figuration produces a larger number of favourable locations for 
maser emission in the corona. 


5. Conclusions 


Using a combination of force-free extrapolation and hydro¬ 
static models, we estimate the co pe /£l e ratio of coronal mag¬ 
netic configurations above active regions. The co pe /£l e ratio is 
important in order to determine in which regime the plasma 
can evolve and what kind of electronic plasma waves can prop¬ 
agate and grow. The main electronic wave modes are the O- 
and X-type for the maser emission, the whistler mode, and the 
upper hybrid mode (electrostatic mode). The different regimes 
in which these modes grow or are sup pres sed is essentially a 
function of the co pe /£l e ratio. Mel rose & Dulkl (1982) have men¬ 



tioned that the growth of the first and second harmonics of the 
X-mode may be responsible for observed radio and hard X-ray 
emissi on in the solar corona. The study of ISharma & Vlahos 
(119841) has defined intervals of co pe /£l e in which the different 
modes are important. For instance, the first harmonic of the X- 
mode has a maximum growth rate for (o pe /£l e < 0.35, whilst the 
first harmonic of the O-mode dominates when 0.35< co pe /£l e 
< 1. The second harmonic of the X-mode is dominant for 1 < 
co p d£l e < 1-45. These values are obtained in the case where 
the coupling with other wave modes is neglected. Therefore, it 
is important to know what values of co pe /£l e can be expected 
in a realistic magnetic field above active regions and thus to 
check if the electron-cyclotron maser emission can be a vi¬ 
able mechanism in the solar corona. To achieve these goals, 
we develop a zero-order magnetohydrostatic model that ne¬ 
glects the coupling between the magnetic field and the plasma 
IRegnier et ali 2008). Applying this method to several active 
regions, we found that co pe /£l e can be less than 0.35 in both 
potential and niff field models; however, this is statistically 
marginal and thus localised in small coronal volumes. 

The two main results obtained from this study are: 


1. The smallest values of co pe /£l e are located where the mag¬ 
netic field strength is the largest at the bottom of the corona 
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Fig. 7. Same as Fig.[T]for AR10486 


over sunspots; this is true whatever the magnetic field 
model used. 

2. Values of (D pe /£l e less than 1 can be found high in the corona 
in the case where highly twisted flux tubes can be found 
within the magnetic configuration; this is only obtained in 
magnetic field models containing electric currents. 


From this new technique for estimating (L> pe /kl e in coronal 
plasmas, we conclude that the maser instability/emission is a 
viable mechanism. As mentioned bv iMelrose & Dulkl (1982), 
the second harmonics of the O- and X-modes are the most pos¬ 
sible modes to be observed, whilst the possibility of observing 
the first harmonic of the X-mode is statistically insignificant. In 
addition, the possible emission is most likely to be localised at 
the bottom of the corona or at coronal heights where the free 
magnetic energy is locally increased (for instance, near twisted 
flux bundles). We also note that both the pressure and magnetic 
field scale-heights are important in order to describe the varia¬ 
tion of the Alfven speed and thus the distribution of aj pe /£l e . 

The low values of oj pe /^e localised high in the corona (see 
Fig. [3] and Fig. [7]) are situated above twisted flux bundles for 


AR8151 and above sheared arcades for AR10846, and are in¬ 
deed related to the complex magnetic topology of the fields 
induced by those structures. To understand better the role that 
these low aj pe /£l e values play in the activity of an active region, 
the next step is to follow the time evolution of an active region 
producing observed radio emission consistent with the maser 
instability. 
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Fig. A.l. Same as Fig.[T]for AR8210 


model, the plasma pressure and density, and the Alfven speed 
are thus given by 


| p(z ) = 

Poex p(-^), 

and 

\ p(z) = 

Poexp(-^p) 

for g = go 

and 



v A (x,y, z) 

B{x,y,z) , 

(Z-Zo\ 

= vR « ext ’ 1 

{ 2 H )’ 


Appendix A: Maser emission and twisted flux 
tubes in a constant gravity field 

For the sake of completeness, we present the E e values de¬ 
rived for AR8151 in the case of a niff field with constant grav¬ 
ity, which exhibits a peculiar feature. For the constant gravity 


where H = IcbT/( pm p go) is the pressure scale-height (kg = 
1.38 10 -23 J K _1 , g = 0.6 for a fully ionised coronal plasma, 
m p = 1.67 10 -27 kg, and go = 274 m s -2 ), and po and po 
are characteristic values of the plasma pressure and density at 
Zo- In Fig. IA.11 an isolated blob of low E e values appears at a 
height of about 130 Mm. Comparing with Fig. 0 the blob is 
not present in the niff model with varying gravity. The blob is 
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located above the highly twisted flux tube that has been defined 


m lRegnier & Amaril d2004l) . A possible explanation for the ap¬ 
pearance of the blob in this particular model is that the exis¬ 
tence of the twisted flux tube changes the decay rate of the mag¬ 
netic field strength (i.e. the magnetic field scale-height), whilst 
the density still follows the hydrostatic decay with height. The 
blob does not appear when the gravity varies with height be¬ 
cause the density drops faster than the magnetic field strength 
as emphasised in Fig. 7 of iRegnier et al.1 ( 2008 ). The blob is 
thus an artefact of this particular model; however, its existence 
emphasises the importance of understanding the variation of 
both the magnetic and pressure scale-heights within the corona. 















